The influence of geometry and topology of quantum graphs on their nonhnear-optical 

properties 
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We analyze the nonlinear optics of quasi one-dimensional quantum graphs and manipulate their 
topology and geometry to generate for the first time nonlinearities in a simple system approaching 
the fundamental limits. Changes in geometry result in smooth variations of the nonlinearities. 
Topological changes between geometrically-similar systems cause profound changes in the nonlinear 
susceptibilities that include a discontinuity due to abrupt changes in the boundary conditions. This 
work may inform the design of new molecules or nano-scale structures for nonlinear optics. 
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I. INTRODUCTION 

A quantum graph is a network of edges and vertices 
upon which particle motion is described by a Hamilto- 
nian with a complete set of eigenstates, an energy spec- 
trum, and flux-conserving boundary conditions An 
electron on the graph is assumed to be tightly bound in 
the transverse direction, yielding a quasi one-dimensional 
dynamical system. When coupled to optical fields, this 
system is a model of nonlinear optical molecular struc- 
tures exhibiting rich topological and geometrical proper- 
ties. 

Quantum graphs with bare and dressed edges have 
been solved using periodic orbit theory and extensively 
studied for their statistical properties and energy spectra 
[l|-0] . This letter is the first to apply quantum graph 
models to nonhnear optics by calculating the hyperpo- 
larizabilities (first and second) using explicit expressions 
for the eigenstates and energy spectra. We employ bare 
graphs which use free-particle states Aij s'm{ksij + 
for each edge connecting vertex pairs j, flux conserva- 
tion at each degree > 1 vertex, and Dirichlet boundary 
conditions at terminal vertices. 

In order to compare the hyperpolarizabilities of quan- 
tum graphs of differing size, we normalize each to the 
fundamental limit [sj to yield the intrinsic first and 
second (7) hyperpolarizability. Unless otherwise stated, 
all values presented here are intrinsic. Secondly, we de- 
fine the molecular x-axis along the direction such that 
Pxxx is maximum. It is typical for ^xxx and "yxxxx to 
be largest along different axes, so the coordinate system 
for describing each may be and usually is different. The 
coordinate axis will be specified only when needed. Pxxx 
ilxxxx) should be understood as the largest x-diagonal 
tensor component of the intrinsic first (second) hyperpo- 
larizability. 

The calculation of the first and second hyperpolariz- 
ability requires a sum over states of products of the dipole 
matrix elements of the graph, divided by energy denom- 



inators ^ . We use a quasi-one dimensional model where 
the electron is tightly confined to move along the edges 
by a transverse well of infinite height and zero width. As 
shown in [loj, the residual effects of the vanishing trans- 
verse dimension are present in the Thomas-Reich-Kuhn 
(TKR) sum rules but drop out of the hyperpolarizabili- 
ties. We verified the correctness of our solutions by cal- 
culating both longitudinal and transverse contributions 
to the sum rules, as well as by showing how they con- 
tribute to the dipole- free sum-over-states formalism , 
so transverse contributions could be omitted from the 
calculation of the hyperpolarizability tensors. 

The first hyperpolarizability f3ijk is a fully-symmetric, 
third rank tensor for planar graphs far from resonance, 
with four nonzero components and an invariant norm 
The second hyperpolarizability is likewise fully- 
symmetric in this limit and has five nonzero components 
and an invariant norm |7|. The norms are given by 
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and contain all of the nonlinear optical physics of the 
model. Each tensor component in Eq. ([T]) may be calcu- 
lated using a sum-over-states expression once the dipole 
moment matrix has been determined, (ol. fl2j 

The computation of the contributions to the dipole mo- 
ment matrix over the edges of the graph is straightfor- 
ward and when summed results in a set of matrix ele- 
ments for the entire graph, from which hyperpolarizabil- 
ities may be calculated with the usual sum-over-states 
expansions, [oi Il2j normalized by their maximum intrin- 
sic values [l^. For the position of the charge in the wire, 
X, we find that 



(2) 
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where i is summed over all edges of the graph, and Oi is 
the angle between the segment i and the external x—axis. 
A similar expression holds for Unm with cos^^ sin^^. 
The factor i^^^ is an integral on edge i of the product of 
the and energy eigenfunctions on the edge times 
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the longitudinal coordinate on the edge, with each i^^^ 
depending solely on the lengths and eigenfunctions on the 
edges. The number of eigenfunctions, their degeneracies 
(if any), and the energy spectra are fixed by the bound- 
ary conditions imposed by the topology of the graph. 
The angular factors in Eq. (j2j) describe the geometry of 
the graph and are identical for graphs with identical ge- 
ometries but different topologies, such as a closed-loop 
triangle graph and a triangle graph with identical shape 
but one open vertex. 

By analyzing classes of graphs with similar geometry 
but different topology, and vice-versa, we can extract 
both the topological and geometrical effects on the non- 
linear optics of the graph. For example, the diagonal 
component of the first hyperpolarizability tensor will al- 
ways take the form 



TABLE I. Effects of geometry and topology in simple graphs. 



Pxxx = ^ (cos 6i COS Oj COS Ojf;) X Aijjf; , 



EjiEr, 
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Equation (|3|) expresses the influence of the angular fac- 
tors describing the geometry of the graph on each seg- 
ment's contribution to the underlying quantum mechan- 
ics of the full graph as embodied in the edge factors i^^^. 
It is thus reasonable to speak of the cosine factors as 
the geometric specifiers. Since the topology of the graph 
determines the boundary conditions on the eigenstates, 
topological effects originate solely in Aijk. 

Aijk has three indices that couple the angular factors in 
ways that are determined by the topology of the graph, 
not its geometry. This means that the differences be- 
tween a closed loop triangle and one with an open vertex 
arises from the quantum states, i.e., the topology of the 
graph. Along the edges of the graph, the electron knows 
nothing about angles; it only knows about the value of 
its wavefunction on any particular edge. 

Table [J displays the four sets of graphs studied. The 
first set are bent wires with different geometries but linear 
topology. The second set are the topologically-equivalent 
closed- loop graphs - again with different geometries. Be- 
low the double line, the next set are graphs with the same 
triangle geometry but which differ in topology. The final 
set are 3-prong irregular (i.e. irrational length ratios) star 
graphs, which are geometrically equivalent, but topolog- 
ically distinct. 

We generated tens of thousands of random configura- 
tions of each type of graph, and calculated the hyperpo- 
larizability tensors for every sample. The simplest graph 
in the first topologically equivalent set is the line segment 
shown in the first row of Table HI and is characterized 
by vanishing Pxxx by virtue of its centrosymmetry, and 
by ^xxxx = —0.126. A bent wire, as shown in the sec- 
ond row of Table [H is topologically equivalent to the line 
but of different geometry (i.e. shape). Figure [1] shows 
the calculated hyperpolarizability for a bent wire that is 
formed by connecting the points (0,0), (x,y) and (1,0). 
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The largest hyperpolarizability is given by = 0.172 
for a wire with ix^y) = (1.85,0). Similarly, the second 
hyperpolarizability varies from -0.126 to +0.007. This 
shows how the hyperpolarizability smoothly varies from 
zero for the straight wire and peaks for a bent wire that 
turns back on itself. For the case of '^xxxxi the best shape 
is given by an acute angle with vertex (0, 0) and endpoints 
(1,-1) and (1.8,0). Despite their simplicity, bent wires 
have an intrinsic first hyperpolarizability that is about 
17% of the fundamental limit. In contrast, most of the 
best molecules fall a factor of 30 short of the limit. [l3[ 

Adding a third segment to a bent wire with no con- 
straint on its orientation yields about the same extremes 
of Pxxx and '^xxxx as the two-wire graph. Imposing geo- 
metrical constraints on the third edge causes minor quan- 
titative change, as can be seen in the line triangles with 
open vertices in Table HI The addition of even more de- 
grees of freedom enables wire configurations of optimum 
shape to be marginally improved and provides limited 
enhancement of nonlinear it ies. For bent wires, ^xxxx can 
be positive or negative. The eigenstates of bent wires are 
non-degenerate, a consequence of the open topology for 
the graph. 

The separation of geometrical and topological fac- 
tors displayed in Eq.(j3j) practically begs for an expan- 
sion of the hyperpolarizabilities into irreducible, spherical 
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FIG. 1. Hyperpolarizability of a bent wire that is defined by 
the coordinates (0,0), (x,y) and (1,0). 



tensors [13, [HI in order to visualize the geometrical and 
topological effects across structures. For fully symmet- 
ric tensors, Pijk is represented by the spherical tensors 
/3i and /3s - corresponding to J = 1 and J = 3 irre- 
ducible representations, while jijki has nonzero 70, 72, 
and 74 - corresponding to J = 0, J = 2, and J = 4 ir- 
reducible representations. Detailed expansions are avail- 
able in the literature [16]. The tensor norms are given by 



33|^ and \-f\^ = |7o|^ + \-f2\^ + \-f4\^; and, 
are identical to those computed in Eq.([T]). 

Consider first the closed-loop topologies from Table 
[D The solutions are a doubly-degenerate set of modes 
with a singlet ground state at zero energy and a peri- 
odic boundary condition on the eigenstates. (Note that 
a complex quadrilateral has the same topology as a sim- 
ple quadrilateral when the crossing edges do not transfer 
probability flux, as we assume here). 

The closed-loop topology severely limits the magni- 
tudes of the first hyperpolarizability across the set of 
geometries. Compared to a simple bent wire, a trian- 
gle loop is a poorer nonlinear optical structure with a 
maximum Pxxx ^ 0.049. More significant, the second 
hyperpolarizability jxxxx is always less than or equal to 
zero in closed-loop graphs. Finally, rapid saturation of 
the nonlinearity with the number of edges occurs, as it 
did for bent wires. Large Pxxx is associated with open, 
isosceles-like shapes, while low-aspect ratio (flat) trian- 
gles yield zero Pxxx- But jxxxx is the most negative for 
flat triangles. Quadrangles and above may be geomet- 
rically squeezed into these shapes, thus explaining how 
the topological features of the loops drive most of the 
physics, while the geometric shapes have only a modest 
effect. 

Figure [2] compares the tensor properties of the distri- 
bution of 10,000 random samples of configurations of the 
closed-loop triangle and quadrangle graphs shown in Ta- 
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FIG. 2. \l3\ vs I3xxx and I7I vs "jxxx for 10,000 randomly- 
sampled graph configurations of fixed topology. Quintangle 
results are not shown, as their tensors look almost identical 
to the quadrangles. 
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FIG. 3. Spherical tensors for \/3\ as a function of the ratio of 
the smallest angle Omin to the largest angle Omax in a closed 
loop (left) and open (right) triangle graph. 



blelH (Note that all plots that follow use 10,000 random 
samples.) The plot on the left shows vs Pxxx^ while 
the plot on the right shows I7I vs jxxxx for the three (fixed 
topology) geometrically-different closed loops. Triangle 
graphs, which have only two angular degrees of freedom, 
have large (absolute) projections onto the x-axis only 
when 1/3 1 is large as well. The quadrilateral, having 3 
degrees of freedom, has the largest value of /3xxx^ but 
marginally so. When has its largest value, Pxxx is not 
optimal. However, is maximum for the same values 
of Pxxx in all geometries. 

The constraint imposed by the closed triangular loop 
on the second hyperpolarizability yields a tight grouping 
of all configurations, as shown by what appears to be 
a smooth red curve. Adding extra degrees of freedom 
yields a greater spread in I7I for a given value of jxxxx^ as 
is also found for the first hyperpolarizability. However, 
"Jxxxx is at its minimum and maximum when I7I is at 
its maximum. When I7I is minimum, "jxxxx appears to 
be minimum only in the limit of an infinite number of 
degrees of freedom. 

In summary, the range of P and 7 over all configura- 
tions is limited by the loop topology, not the geometry. 

Figure [3] shows a plot of the contributions to \/3\ from 
the two spherical tensors for a closed triangle and an open 
vertex triangle. For closed loops, large graphs are 
dominated by the tensor when Omin = Omax implying 
that /3 is maximal for an equilateral triangle. The black 
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FIG. 4. Topological dependence for fixed geometry. The fig- 
ures show the variation in the hyperpolarizability tensor com- 
ponents and their norms for open and closed triangle graphs. 



points are for isosceles triangles. The conclusion is that 
high aspect ratio triangles have near-zero /3xxx regardless 
of their orientation. 

For triangles with an open vertex, corresponding to 
a change in topology from the closed triangle, the re- 
lationship between the spherical tensor components and 
the angle ratio are profoundly different. /3i exceeds /3s 
and 1/3 1 is substantially larger for the open vertex trian- 
gles. Furthermore, isosceles triangles no longer give the 
largest P and the equilateral triangle has smallest 
compared with largest in the case of the closed trian- 
gle. Topology therefore has a profound effect on both the 
magnitude and character of the first hyperpolarizability. 

The profound effect of topology on the first and sec- 
ond hyperpolarizability is emphasized in Figure |H For 
the closed loop, both /3xxx and \/3\ are smaller; and, the 
range of (3xxx is substantially less than for When 
a vertex is opened, /3xxx and \f3\ are both larger, and 
there exists a configuration for which pxxx = \P\ for all 
The opposite is true for 7. The closed configuration 
yields the largest value of I7I and jxxxx- Additionally, the 
same span of geometrical configurations for each leads to 
divergent behavior as shown in the figures. 

Pi is always smaller than /^s in the closed triangle but 
are both comparable in the open one. The second hyper- 
polarizability 7 is dominated by 72 in all cases, and each 
ji in the loop topology approximately parallels 7^ in the 
open vertex case; but, the values of 7^ are larger for the 
loops. 

These observations can be understood as follows. 
Opening a single vertex in a closed triangle removes half 
the eigenstates and shifts the ground state to nonzero 
energy. The result is a shift in the vector component 
of 1/3 1 from an insignificant contributor in closed trian- 
gles to a substantial contributor for the open triangles. 
The topological change causes two significant changes in 



Jxxxx^ as it allows for configurations with exactly zero 
norm as well as shapes with positive jxxxx for open tri- 
angles. Though geometrically identical, the closed and 
open topologies differ by over a factor of three in /3xxx^ 
and the open triangles have configurations with positive 

'Jxxxx- 

For simple quadrangles, the distribution of 7 fills a 
greater area in the \f3j\ — \/3\ plane than for triangles, 
reflecting the additional degree of freedom in its con- 
figuration space. In a triangle, two angles define a set 
of similar triangles - all having the same geometry and 
same intrinsic nonlinear it ies. In a quadrangle, three an- 
gles do not form a unique similarity class, leading to a 
greater spread of intrinsic nonlinearity. The additional 
degrees of freedom afforded to a general quadrangle fills 
even a greater part of the plane, yet /3xxx increases only 
modestly with increasing degrees of freedom. The loop 
topology constrains Pxxx to values well under that of the 
simple bent wires and keeps jxxxx negative for any shape. 

Table [J summarizes these results and compares the 
closed-loop topology with three line configurations, each 
of which is a bent wire or a composite of bent wires dis- 
cussed previously. Opening a single vertex in a closed- 
loop triangle eliminates the eigenstate degeneracies and 
produces a ground-state with nonzero energy. The bent 
wire triangle with one open vertex and no flux circula- 
tion has a much larger Pxxx 0.133, and a positive, 
maximum jxxxx- In these triangle graphs, the extremes 
of Pxxx and jxxxx occur for the opposite geometries of 
those in the closed-loop triangle. 

To quantify the geometrical effect of the 'openness' of 
a loop, we define the dimensionless geometrical parame- 
ter ^ as the ratio of the area of a loop to the square of its 
perimeter. Since a polygon with all edge lengths equal to 
each other and all angles equal to each other has maxi- 
mum ^, we normalize all ^ to this value to get the intrinsic 
geometric parameter. Fig. [5] shows \f3\ and I7I (normal- 
ized to unity) as a function of ^ for both closed (loop) 
and open (bent wire, one vertex open) triangles. The 
black points in the Figure are isosceles-like triangles and 
^ = 1 is an equilateral triangle. Opening a second vertex 
in the bent-wire triangle increases the Pxxx toward that 
of a bent wire while still yielding positive jxxxx configu- 
rations. Opening all vertices produces three single wires, 
for which Pxxx is exactly zero and jxxxx approaches the 
value it has for a single bent wire. In such systems, the 
eigenstates are constructed from three 1-wire graphs by 
ordering the individual energies of each wire to form the 
composite and demanding that superpositions maintain 
unitarity. 

The final set of graphs in Table [H the so-called star 
graphs, are fundamentally different than any discussed 
so far. The 3-prong star graph is an irregular graph with 
four degrees of freedom (two angles and the lengths of two 
edges), providing much larger nonlinearities than any of 
the regular graphs. The spectra of the 3-prong topology 
has been solved in closed form using a periodic orbit ex- 
pansion [13, and properties of the eigenfunctions have 
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FIG. 5. Scaling of (left) and I7I (right), both normalized to 
unity, as a function of geometric parameter ^, for closed- and 
open-loop topologies. The black points represent isosceles-like 
geometries. 




FIG. 6. Topological dependence for fixed geometry of star 
graphs, showing on the left (right) a very large enhancement 

in I^XXX i^'^XXXX 

) for the closed star with three attached prongs 
(red) compared to the same graph with one prong detached 
(green) . 

been studied [18]. For nonlinear optics, a numerical solu- 
tion of the secular equation for the star graph wave equa- 
tion provides the simplest way to get the energies and 
eigenfunctions required to compute the dipole moment 
matrix for use in the sum-over-states construction of Pxxx 



and ^xxxx- The result is a maximum Pxxx of order 0.55, 
over half of the fundamental limit and the largest calcu- 
lated to date for any simple structure. Fig [6] illustrates 
the dramatic change caused by a simple topological shift 
in a graph from three connected prongs to two. To the 
extent that stars are effective nano-wire[19] and molec- 
ular models, they cou ld p rovide the molecular designer 
and nano-technologist |2Qj with guidelines for achieving 
breakthroughs in nonlinear optics. As expected, a ma- 
jor change in the topology of the star graph by opening 
a prong at the central vertex changes the highly-active 
3-fork into a graph resembling a bent wire system. An 
obvious conclusion of this analysis is that combinations 
of simple structures such as bent wires and stars may lead 
to structures with even larger nonlinear optical response. 
More complex structures are under study. 

In summary, we have presented for the first time an 
exact, quantitative dependence of the nonlinear optics of 
quantum graphs on their geometry and topology. We 
have shown that the effects of the topology of geomet- 
rically similar graphs dominate those of the geometry 
of topologically similar graphs. Topology largely deter- 
mines the eigenstates and spectra, whereas the geometry 
mainly affects the projections of the graph edges in the 
x—y plane. Closed loop graphs always have non-optimum 
I3xxx and negative ^xxxx-, but opening a vertex immedi- 
ately raises Pxxx by over a factor of three and enables 
graphs with positive '^xxxx- We have also verified that 
additional degrees of freedom enhance the nonlinearity of 
the graph, unless a fundamental topological constraint is 
in place, such as a closed loop boundary condition on the 
eigenstates. Finally, we have shown that the star graph 
has the largest intrinsic hyperpolarizability calculated to 
date in a simple quantum system. 

SS, JS and MGK thank the National Science Foun- 
dation (ECCS-1 128076) for generously supporting this 
work. 
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